Transcriptomic analysis of skin immunity genes in the Chinese spiny frog (Quasipaa spinosa) after Proteus mirabilis infection.

Recently, populations of Chinese spiny frogs (Quasipaa spinosa), an important amphibian species in China, have decreased, mainly due to a disease caused by the gram-negative bacteria Proteus mirabilis. To elucidate the immune response of the frogs, this study aimed to identify novel candidate genes functionally associated with P. mirabilis infection-induced "rotting skin" disease. Chinese spiny frogs were infected with P. mirabilis, and the skin transcriptome was sequenced using the MGISEQ-2000 platform. A total of 233,965 unigenes were obtained by sequencing, of which 27.23 % were known genes. Screening of differentially expressed genes (DEGs) indicated 210 unigenes differentially expressed after P. mirabilis infection, of which 132 unigenes were up-regulated, and 78 unigenes were down-regulated. Using Kyoto Encyclopedia of Genes and Genomes enrichment analysis, DEGs were identified as enriched in signal pathways, such as oxidative phosphorylation, apoptosis, and the Janus kinase-signal transducer and activator of transcription pathway. Of the DEGs, there was a significant upregulation of the colony stimulating factor 2 receptor beta common subunit, interleukin 2 receptor subunit gamma, cathelicidin antimicrobial peptide, interleukin-17 receptor E, receptor-interacting serine/threonine-protein kinase 3, and pulmonary surfactant-associated protein D immune genes following P. mirabilis infection. Conversely, scavenger receptor cysteine-rich domain-containing group B protein, tumor protein p53 inducible nuclear protein 2, suppressor of cytokine signaling 2, and metalloreductase STEAP3 were significantly downregulated. In conclusion, the first skin transcriptome database of Chinese spiny frogs was established, and several immune genes were identified to elucidate the pathogenic mechanism of "skin rot" in Chinese spiny frogs and other cultured frogs.

Karyotypes of water frogs from the Pelophylax esculentus complex: results of cross-species chromosomal painting.

Amphibian species have the largest genome size enriched with repetitive sequences and relatively similar karyotypes. Moreover, many amphibian species frequently hybridize causing nuclear and mitochondrial genome introgressions. In addition, hybridization in some amphibian species may lead to clonality and polyploidization. All such events were found in water frogs from the genus Pelophylax. Among the species within the genus Pelophylax, P. esculentus complex is the most widely distributed and well-studied. This complex includes two parental species, P. ridibundus and P. lessonae, and their hybrids, P. esculentus, reproducing hemiclonally. Parental species and their hybrids have similar but slightly polymorphic karyotypes, so their precise identification is still required. Here, we have developed a complete set of 13 chromosome painting probes for two parental species allowing the precise identification of all chromosomes. Applying chromosomal painting, we identified homologous chromosomes in both parental species and orthologous chromosomes in their diploid hemiclonal hybrids. Comparative painting did not reveal interchromosomal exchanges between the studied water frog species and their hybrids. Using cross-specific chromosome painting, we detected unequal distribution of the signals along chromosomes suggesting the presence of species-specific tandem repeats. Application of chromosomal paints to the karyotypes of hybrids revealed differences in the intensity of staining for P. ridibundus and P. lessonae chromosomes. Thus, both parental genomes have a divergence in unique sequences. Obtained chromosome probes may serve as a powerful tool to unravel chromosomal evolution in phylogenetically related species, identify individual chromosomes in different cell types, and investigate the elimination of chromosomes in hybrid water frogs.

Expression of cholesterol synthesis and steroidogenic markers in females of the Chinese brown frog (Rana dybowskii) during prespawning and prehibernation.

The oviduct of the Chinese brown frog (Rana dybowskii) expands in prehibernation rather than in prespawning, which is one of the physiological phenomena that occur in the preparation for hibernation. Steroid hormones are known to regulate oviductal development. Cholesterol synthesis and steroidogenesis may play an important role in the expansion of the oviduct before hibernation. In this study, we investigated the expression patterns of the markers that are involved in the de novo steroid synthesis pathway in the oviduct of R. dybowskii during prespawning and prehibernation. According to histological analysis, the oviduct of R. dybowskii contains epithelial cells, glandular cells, and tubule lumens. During prehibernation, oviductal pipe diameter and weight were significantly larger than during prespawning. 3-Hydroxy-3-methylglutaryl CoA reductase (HMGCR), low-density lipoprotein receptor (LDLR), steroidogenic acute regulatory protein (StAR), cytochrome P450 cholesterol side-chain cleavage enzyme (P450scc), and steroidogenic factor 1 (SF-1) were detected in epithelial cells in prehibernation and glandular cells during prespawning. HMGCR, LDLR, StAR, and P450scc protein expression levels were higher in prehibernation than during prespawning, but the SF-1 protein expression level did not significantly differ. HMGCR, LDLR, StAR, P450scc (CYP11A1), and SF-1 (NR5A1) mRNA expression levels were significantly higher in prehibernation compared with prespawning. The transcriptome results showed that the steroid synthesis pathway was highly expressed during prehibernation. Existing results indicate that the oviduct is able to synthesize steroid hormones using cholesterol, and that steroid hormones may affect the oviductal functions of R. dybowskii.

Subchronic toxicity of iron-selenium nanoparticles on oxidative stress response, histopathological, and nuclear damage in amphibian larvae Rana saharica.

In this work, we evaluated the subchronic toxicity of FeSe nanoparticles (NPs) in tadpoles of Rana saharica. Tadpoles were exposed for 1-3 weeks to FeSe NPs at 5 mg/L and 100 mg/L rates. Parameters of oxidative stress were measured in whole larvae, and the micronucleus test was performed on circulating blood erythrocytes. We noted a disturbance of the detoxification systems. Enzymatic and non-enzymatic data showed that exposure to FeSe NPs involved a highly significant depletion of GSH, a significant increase in GST activity, and a lipid peroxidation associated with a highly significant increase in MDA. We also noted a neurotoxic effect characterized by a significant inhibition of AChE activity. A micronucleus test showed concentration-dependent DNA damage. This research reveals that these trace elements, in their nanoform, can cause significant neurotoxicity, histopathologic degeneration, cellular and metabolic activity, and genotoxic consequences in Rana larvae.

EThcD Benefits for the Sequencing Inside Intramolecular Disulfide Cycles of Amphibian Intact Peptides.

Disulfide bonds formed by a pair of cysteine residues in the peptides' backbone represent a certain problem for their sequencing by means of mass spectrometry. As a rule, in proteomics, disulfide bonds should be cleaved before the analysis followed by some sort of chemical derivatization. That step is time-consuming and may lead to losses of minor peptides of the analyzed mixtures due to incomplete reaction, adsorption on the walls of the vials, etc. Certain problems in the de novo top-down sequencing of amphibian skin peptides are caused by the C-terminal disulfide loop, called the Rana box. Its reduction with or without subsequent derivatization was considered to be an unavoidable step before mass spectrometry. In the present study, EThcD demonstrated its efficiency in sequencing intact disulfide-containing peptides without any preliminary derivatization. Applied to the secretion of three frog species, EThcD provided the full sequence inside the intramolecular disulfide cycle for all S-S-containing peptides found in the samples, with the only exception being diarginine species. Proteolytic fragments, which are shorter than the original peptides, were helpful in some cases. HCD should be mentioned as a complementary tool to the EThcD tool, being useful as a confirmation method for some sequence details.

Decreased acetylation of HDGF improves oviduct production in Rana dybowskii, Rana amurensis, and Rana huanrenensis.

The oviduct of female Rana dybowskii is a functional food and can be used as a component of Traditional Chinese medicine. The differentially expressed genes enriched was screened in cell growth of three Rana species. We quantitatively analyzed 4549 proteins using proteomic techniques, enriching the differentially expressed proteins of Rana for growth and signal transduction. The results showed that log2 expression of hepatoma-derived growth factor (HDGF) was increased. We further verified 5 specific differential genes (EIF4a, EIF4g, HDGF1, HDGF2 and SF1) and found that HDGF expression was increased in Rana dybowskii. Through acetylation modification analysis, we identified 1534 acetylation modification sites in 603 proteins, including HDGF, and found that HDGF acetylation expression was significantly reduced in Rana dybowskii. Our results suggest that HDGF is involved in the development of oviductus ranae, which is regulated by acetylation modification.

Hepatic citrate synthase suppression in the freeze-tolerant wood frog (Rana sylvatica).

The wood frog, Rana sylvatica endures whole body freezing for weeks/months while overwintering at subzero temperatures. Survival of long-term freezing requires not only cryoprotectants but also strong metabolic rate depression (MRD) and reorganization of essential processes in order to maintain a balance between ATP-producing and ATP-consuming processes. Citrate synthase (CS) (E.C. 2.3.3.1) is an important irreversible enzyme of the tricarboxylic acid (TCA) cycle and forms a crucial checkpoint for many metabolic processes. Present study investigated the regulation of CS from wood frog liver during freezing. CS was purified to homogeneity by a two-step chromatographic process. Kinetic and regulatory parameters of the enzyme were investigated and, notably, demonstrated a significant decrease in the Vmax of the purified form of CS from frozen frogs as compared to controls when assayed at both 22 °C and 5 °C. This was further supported by a decrease in the maximum activity of CS from liver of frozen frogs. Immunoblotting also showed changes in posttranslational modifications with a significant decrease in threonine phosphorylation (by 49 %) for CS from frozen frogs. Taken together, these results suggest that CS is suppressed and TCA flux is inhibited during freezing, likely to support MRD survival of harsh winters.

Effects of cadmium and diethylhexyl phthalate on skin microbiota of Rana chinensis tadpoles.

Skin microbiotas play a crucial role in the health, homeostasis, and immune function of amphibians. The contaminants in water could affect the structure and composition of microbial communities. The effects of coexisting pollutants on frogs cannot be adequately explained by a single exposure due to the coexistence of Cd and DEHP in the environment. Following exposure to Cd and/or DEHP, we examined the histological characteristics of Rana chensinensis tadpoles. We also used the 16S rRNA gene sequencing technique to assess the relative abundance of skin microbial communities among tadpoles from each treatment group. Our findings indicate that R. chensinensis' skin experienced some degree of injury due to exposure to Cd and DEHP, which led to the imbalance of their skin microbial community homeostasis and thus interfered with the normal trial status of the host. That may eventually lead to the decline of the amphibian population.

Temporins: Multifunctional Peptides from Frog Skin.

Temporins are short peptides secreted by frogs from all over the world. They exert antimicrobial activity, mainly against Gram-positive bacteria, including resistant pathogens; recent studies highlight other possible applications of these peptides as anticancer or antiviral agents. This review is meant to describe the main features of temporins produced by different ranid genera. Due to the abundance of published papers, we focus on the most widely investigated peptides. We report studies on their mechanism of action and three-dimensional structure in model systems mimicking bacterial membranes or in the presence of cells. The design and the antimicrobial activity of peptide analogues is also described, with the aim of highlighting elements that are crucial to improve the bioactivity of peptides while reducing their toxicity. Finally, a short section is dedicated to the studies aimed at applying these peptides as drugs, to produce new antimicrobial materials or in other technological uses.

Diversification dynamics in Caribbean rain frogs (Eleutherodactylus) are uncoupled from the anuran community and consistent with adaptive radiation.

Adaptive radiation is proposed to play a key role in generating differences in species richness among lineages and geographical regions. Due to the importance of ecological divergence in adaptive radiation, species richness is predicted to be influenced by equilibrium diversity dynamics, although the concept continues to generate much debate. An additional important question is whether radiating clades have intrinsic biological characteristics that make them particularly prone to diversify. We tackle these questions by analysing (i) the temporal patterns of diversification of Caribbean Eleutherodactylus frogs, and (ii) assembly of the complete native anuran community of the Caribbean archipelago (197 species), testing for the presence of equilibrium dynamics and whether diversification patterns of Eleutherodactylus differ from those of the rest of the Caribbean anurans. Diversification rates follow the predicted pattern of rapid diversification early in the radiation which gradually decreases towards the present. Eleutherodactylus diversification is significantly faster than that of the Caribbean anuran community, and although equilibrium dynamics influence richness of all Caribbean anurans, Eleutherodactylus shows higher carrying capacity. Our results indicate that ecological opportunity per se is not sufficient for adaptive radiation and that diverse lineages present intrinsic characteristics that enable them to make the most of available opportunity.

Evaluation of the Synthetic Multifunctional Peptide Hp-MAP3 Derivative of Temporin-PTa.

In recent years, antimicrobial peptides isolated from amphibian toxins have gained attention as new multifunctional drugs interacting with different molecular targets. We aimed to rationally design a new peptide from temporin-PTa. Hp-MAP3 (NH2-LLKKVLALLKKVL-COOH), net charge (+4), hydrophobicity (0.69), the content of hydrophobic residues (69%), and hydrophobic moment (0.73). For the construction of the analog peptide, the physicochemical characteristics were reorganized into hydrophilic and hydrophobic residues with the addition of lysines and leucines. The minimum inhibitory concentration was 2.7 to 43 µM against the growth of Gram-negative and positive bacteria, and the potential for biofilm eradication was 173.2 µM. Within 20 min, the peptide Hp-MAP3 (10.8 µM) prompted 100% of the damage to E. coli cells. At 43.3 µM, eliminated 100% of S. aureus within 5 min. The effects against yeast species of the Candida genus ranged from 5.4 to 86.6 µM. Hp-MAP3 presents cytotoxic activity against tumor HeLa at a concentration of 21.6 µM with an IC50 of 10.4 µM. Furthermore, the peptide showed hemolytic activity against murine erythrocytes. Structural studies carried out by circular dichroism showed that Hp-MAP3, while in the presence of 50% trifluoroethanol or SDS, an α-helix secondary structure. Finally, Amphipathic Hp-MAP3 building an important model for the design of new multifunctional molecules.

Effects of minimum Strigea robusta (Digenea: Strigeidae) cercariae doses and localization of cysts on the anomaly P manifestation in Pelophylax lessonae (Anura: Ranidae) tadpoles.

The anomaly P is a mass morphological anomaly found in some populations of anuran amphibians (water frogs of the genus Pelophylax and toads of the genera Bufo and Bufotes) caused by the parasitic flatworm Strigea robusta. Minimum dose of cercariae for the appearance of the anomaly P remains unknown. However, it is important information for understanding of host population dynamics after invasion and the effects of the parasite on the second intermediate hosts. Herein, the invasion properties of S. robusta in Pelophylax lessonae tadpoles (Anura: Ranidae) and minimum dose for appearance of mild and severe forms of the anomaly P syndrome were described after direct experiments with certain numbers of cercariae exposure. Experimental groups of tadpoles have been exposed to eight doses of cercariae (2, 4, 6, 8, 10, 12, 14 and 16). A total of 63.8% tadpoles survived to the end of this experiment. It was revealed that a mild form of the anomaly P (polydactyly) can appear after infection by two cercariae, while the severe form traits appear after infection by four cercariae. The mean number of detected encysted metacercariae was reached to 53.5%. Differences in infection rates can be explained by the presence of an individual immune response in tadpoles or by the presence of different genetic lineages of the parasite infecting the same snail, which have different infectious potential. Low doses of infection leading to the induction of anomalies characterize S. robusta as a highly pathogenic species for amphibian species that are susceptible to infection and show an abnormal phenotype.

mRNA m6 A methylation in wood frog brain is maintained during freezing and anoxia.

Freeze tolerance is an adaptive strategy that wood frogs (Rana sylvatica) use to survive the subzero temperatures of winter. It is characterized by a variety of metabolic and physiological changes that facilitate successful freezing and anoxia. As both mRNA regulation and posttranslation protein modification have been implicated in freeze tolerance, we hypothesized that posttranslational RNA regulation is also involved in coordinating freeze-thaw cycles and metabolic rate depression. As such, we investigated the most abundant RNA modification, adenosine methylation (N6 -methyladenosine; m6 A) in wood frog brains during 24 h periods of freezing and anoxia. This was followed by an examination of levels of RNA methyltransferases, demethyltransferases, and the readers of RNA methylation. Despite relative levels of methylation on mRNA remaining constant throughout freezing and anoxia, a significant increase in relative abundance of m6 A methyltransferases METTL3 and METTL14 was observed. In addition, we investigated the effect of m6 A RNA methylation on mRNA triaging to stress granules and report a significant increase in stress granule markers TIAR and TIA-1 in both freezing and anoxia. Our findings are the first report of RNA posttranslational regulation during metabolic rate depression in the wood frog brain and suggest that the dynamic RNA methylation observed is not directly linked to mRNA regulation during periods of extreme metabolic reorganization, warranting future investigations.

Freeze-induced suppression of pyruvate kinase in liver of the wood frog (Rana sylvatica).

The wood frog (Rana sylvatica) undergoes physiological and metabolic changes to withstand subzero temperatures and whole body freezing during the winter months. Along with metabolic rate depression, high concentrations of glucose are produced as a cryoprotectant by liver and distributed to all other tissues. Pyruvate kinase (PK; EC:2.7.1.40), the final enzyme of glycolysis, plays an important role in the modulation of glucose metabolism and, therefore, overall metabolic regulation. The present study investigated the functional and kinetic properties of purified PK from liver of control (5 °C acclimated) and frozen (-2.5 °C for 24 h) wood frogs. Liver PK was purified to homogeneity by a two-step chromatographic process, followed by analysis of enzyme properties. A significant decrease in the affinity of PK for its substrates, phosphoenolpyruvate (PEP) and adenosine diphosphate (ADP) at 22 °C and 5 °C was noted in liver from frozen frogs, as compared with controls. Immunoblotting also revealed freeze-responsive changes in posttranslational modifications with a significant increase in serine and threonine phosphorylation by 1.46-fold and 1.73- fold for PK from frozen frogs as compared with controls. Furthermore, a test of thermal stability showed that PK from liver of frozen wood frogs showed greater stability as compared with PK from control animals. Taken together, these results suggest that PK is negatively regulated, and glycolysis is suppressed, during freezing. This response acts as an important survival strategy for maintaining continuously elevated levels of cryoprotectant in frogs while they remain in a hypometabolic frozen state.

Perfluorooctanoic acid and perfluorooctanesulfonic acid induce immunotoxicity through the NF-κB pathway in black-spotted frog (Rana nigromaculata).

Perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS) are widely detected in the environment and wild animals, thus posing a threat to wildlife and public health; however, knowledge about their immunotoxicity and the underlying mechanism remains limited. In the present study, male black-spotted frogs (Rana nigromaculata) were exposed to environmentally relevant concentrations (0, 1, and 10 µg/L) of PFOA or PFOS for 21 days; subsequently, biochemical analysis, molecular docking, and gene expression determination were conducted. The results indicated that exposure to 10 µg/L PFOA decreased the serum levels of immunoglobulin A. PFOS exposure significantly increased the hepatic levels of interleukin-1ß, interleukin-6, tumor necrosis factor-α, interferon-γ, and nitric oxide; but PFOA significantly increased the levels of only tumor necrosis factor-α. Furthermore, PFOA and PFOS exposure significantly decreased the activity of inducible nitric oxide synthase and total nitric oxide synthase. IBRv2 analysis indicated that PFOA and PFOS had a similar effect on these immune indicators, but PFOS was more toxic than PFOA. Molecular docking revealed that PFOA and PFOS can bind to nuclear factor-κB (NF-κB) by forming stable hydrogen bonds. PFOA and PFOS exposure upregulated the gene expression of NF-κB and its downstream genes. Significant correlations between the expression of genes involved in the NF-κB pathway and immune-related indicators suggests that PFOA- and PFOS-induced immunotoxicity was associated with the activation of NF-κB. Our findings provide novel insights into the potential role of NF-κB in immunotoxicity induced by PFOA and PFOS in frogs.

The regulation of m6A-related proteins during whole-body freezing of the freeze-tolerant wood frog.

Rana sylvatica (also known as Boreorana sylvatica) is one of the few vertebrates that spend extreme winters showing no physiological signs of life. Up to 70% of the total body water of the wood frog freezes as extracellular ice. Survival in extreme conditions requires regulation at transcriptional and translational levels to activate prosurvival pathways. N6-methyladenosine (m6A) methylation is one of the most common RNA modifications, regulating transcript processing and translation by executing important functions that affect regulatory pathways in stress conditions. In the study, regulation of m6A-related proteins in the liver of R. sylvatica was analyzed during 24 h frozen and 8 h thaw conditions. Decreases in the activity of demethylases of 28.44 ± 0.4% and 24.1 ± 0.9% of control values in frozen and thaw tissues, respectively, were observed. Total protein levels of m6A methyltransferase complex components methyltransferase-like 14 and Wilm's tumor associated protein were increased by 1.28-fold and 1.42-fold, respectively, during freezing. Demethylase fat mass and obesity, however, showed a decreasing trend, with a significant decrease in abundance during recovery from frozen conditions. Levels of mRNA degraders YTHDF2 and YTHDC2 also decreased under stress. Overall, increased levels of m6A methylation complex components, and suppressed levels of readers/erasers, provide evidence for the potential role of RNA methylation in freezing survival and its regulation in a hypometabolic state.

Comparison of transcriptome responses of the liver, tail fin, and olfactory epithelium of Rana [Lithobates] catesbeiana tadpoles disrupted by thyroid hormones and estrogen.

Thyroid hormones (THs) are important developmental regulators in vertebrates, including during the metamorphosis of a tadpole into a frog. Metamorphosis is a post-embryonic developmental period initiated by TH production in the tadpole thyroid gland. The two main bioactive forms of TH are L-thyroxine (T4) and 3,5,3'-triiodothyronine (T3); these hormones have overlapping but distinct mechanisms of action. Premetamorphic tadpoles are highly responsive to TH and can be induced to metamorphose through exogenous TH exposure, making them an important model for both the study of vertebrate TH signaling and endocrine disrupting chemicals (EDCs). It is important to differentiate TH-mediated responses from estrogenic responses in premetamorphic tadpoles when assessing dysregulation by EDCs as crosstalk between the two endocrine systems is well-documented. Herein, we compare the RNA-sequencing-derived transcriptomic profiles of three TH-responsive tissues (liver, olfactory epithelium, and tail fin) in premetamorphic bullfrog (Rana [Lithobates] catesbeiana) tadpoles exposed to T3, T4, and estradiol (E2). These profiles were generated using the latest available genome assembly for the species. The data indicate that there is a clear distinction, and little overlap, between the transcriptomic responses elicited by E2 and the THs. In contrast, within the THs, the T3- and T4-induced transcriptomic profiles generally show considerable overlap; however, the degree of overlap is highly tissue-dependent, illustrating the importance of distinguishing the two THs and the affected signaling pathways within the target tissue type when evaluating hormone active agents. The data herein also show that E2 and TH treatment can uniquely induce significant changes in expression of their respective "classic" bioindicator transcripts vtg (E2) and thra, thrb, and thibz (THs). However, care must be taken in the interpretation of increased vep or esr1 transcripts as a change in transcript levels can be induced by THs rather than solely E2.

First Evidence of Anti-Steatotic Action of Macrotympanain A1, an Amphibian Skin Peptide from Odorrana macrotympana.

Many different amphibian skin peptides have been characterized and proven to exert various biological actions, such as wound-healing, immunomodulatory, anti-oxidant, anti-inflammatory and anti-diabetic effects. In this work, the possible anti-steatotic effect of macrotympanain A1 (MA1) (FLPGLECVW), a skin peptide isolated from the Chinese odorous frog Odorrana macrotympana, was investigated. We used a well-established in vitro model of hepatic steatosis, consisting of lipid-loaded rat hepatoma FaO cells. In this model, a 24 h treatment with 10 µg/mL MA1 exerted a significant anti-steatotic action, being able to reduce intracellular triglyceride content. Accordingly, the number and diameter of cytosolic lipid droplets (LDs) were reduced by peptide treatment. The expression of key genes of hepatic lipid metabolism, such as PPARs and PLINs, was measured by real-time qPCR. MA1 counteracted the fatty acid-induced upregulation of PPARγ expression and increased PLIN3 expression, suggesting a role in promoting lipophagy. The present data demonstrate for the first time a direct anti-steatotic effect of a peptide from amphibian skin secretion and pave the way to further studies on the use of amphibian peptides for beneficial actions against metabolic diseases.

Assessing the hepatotoxicity of PFOA, PFOS, and 6:2 Cl-PFESA in black-spotted frogs (Rana nigromaculata) and elucidating potential association with gut microbiota.

Pollution caused by per- and polyfluoroalkyl substances (PFASs) has become a major global concern. The association between PFAS-induced hepatotoxicity and gut microbiota in amphibians, particularly at environmentally relevant concentrations, remains elusive. Herein we exposed male black-spotted frogs (Rana nigromaculata) to 1 and 10 µg/L waterborne perfluorooctanoic acid (PFOA), perfluorooctanesulfonic acid (PFOS), and 6:2 chlorinated polyfluorinated ether sulfonate (6:2 Cl-PFESA) for 21 days; subsequently, liver histopathological, oxidative stress, molecular docking, gene/protein expression, and gut microbiome analyses were conducted. PFOS and 6:2 Cl-PFESA exposure enhanced serum alanine aminotransferase and aspartate aminotransferase activities, and markedly increased hepatic area of vacuoles and inflammatory cell infiltration, while PFOA exposure increased serum alanine aminotransferase but not aspartate aminotransferase activities and affected hepatic area of vacuoles and inflammatory cell infiltration to a lesser extent. All three PFASs elevated catalase, glutathione S-transferase, and glutathione peroxidase activities and glutathione and malondialdehyde contents in the liver, suggesting the induction of oxidative stress. Further, PFASs could bind to mitogen-activated protein kinases (p38, ERK, and JNK), upregulating not only their expression but also the expression of downstream oxidative stress-related genes and that of P-p38, P-ERK, and Nrf2 proteins. In addition, PFAS exposure significantly increased the relative abundance of Proteobacteria and Delftia and decreased that of Firmicutes and Dietzia, Mycoplasma, and Methylobacterium-Methylorubrum in the order of PFOS ≈ 6:2 Cl-PFESA > PFOA. Altogether, it appears that PFOS and 6:2 Cl-PFESA are more toxic than PFOA. Finally, microbiota function prediction, microbiota co-occurrence network, and correlation analysis between gut microbiota and liver indices suggested that PFAS-induced hepatotoxicity was associated with gut microbiota dysbiosis. Our data provide new insights into the role of gut microbiota in PFAS-induced hepatotoxicity in frogs.